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RAPID  ELECTRON  BEAM  ACCELERATORS 
(REBA-TRONS) 

I.  Introduction 

Ultra-high  current  accelerators  are  rapidly  becoming  an  active  area 
I  2 

of  research.  ’  The  development  of  these  devices  is  mainly  motivated  by 
a  variety  of  potential  applications1 *J *  that  are  extended  over  several 
areas,  including  environment,  food  processing,  radiation  sources,  x-ray 
radiography  and  national  defense. 

Among  the  various  accelerating  schemes  that  have  the  potential  to 
produce  ultra-high  power  electron  beams,  induction  accelerators^ *^ 
appear  to  be  the  most  promising.  Induction  accelerators  are  inherently 
low  impedance  devices  and  thus  are  ideally  suited  to  drive  high  current 
beams.  The  acceleration  process  is  based  on  the  inductive  electric 
field  produced  by  a  time  varying  magnetic  field. 

5-13 

Quite  naturally,  induction  accelerators  are  divided  into  linear 
14-19 

and  cyclic  devices.  In  linear  devices  the  accelerating  field  is 

localized  in  the  gap,  while  in  their  cyclic  counterparts  the  electric 

field  is  continuous  along  the  orbit  of  the  accelerated  particles.  Both 

cyclic  and  linear  devices  require  the  same  total  magnetic  flux  change  to 

achieve  a  given  energy  increment.  However,  in  linear  accelerators  the 

total  change  of  flux  occurs  in  one  transit  time,  typically  in  less  than 

100  nsec,  while  in  cyclic  accelerators  the  same  change  occurs  over 

several  thousand  revolutions  in  a  typical  time  of  one  msec. 

As  a  consequence  of  the  slow  acceleration,  the  accelerated  beam 

must  be  confined  by  the  focusing  magnetic  field  over  long  periods  of 

time  and  thus  field  errors,  instabilities  and  radiation  losses  impose 

limitations  on  the  cyclic  accelerators.  These  limitations  can  be 

substantially  relaxed  if  the  acceleration  could  occur  rapidly  as  in 
Manuscript  approved  February  7,  1985. 


linear  accelerators.  Therefore,  a  device  that  combines  the  rapid 
acceleration  of  linear  accelerators  and  the  compact  size  of  cyclic 

accelerators  is  highly  desirable.  In  this  paper,  we  propose  such  a 
hybrid  scheme  that  combines  most  of  the  advantages  of  linear  and  cyclic 
accelerators.  This  device  has  been  named  REBA-TRON  (Rapid  Electron  Beam 
Accelerator).  The  rebatron  is  similar  to  the  racetrack  Induction 
accelerator.  u 

A  rebatron  is  shown  schematically  in  Fig.  1.  The  high  gradient 

localized  field  that  is  responsible  for  the  rapid  acceleration  is 

produced  by  convoluted  parallel  transmission  lines,  although,  other 

transmission  lines  may  be  more  appropriate  in  an  actual  system.  Since 

the  acceleration  occurs  over  a  few  usee,  the  constraints  imposed  on  the 

vertical  field  are  very  stringent.  In  an  actual  device  the  vertical 

field  is  generated  by  two  coaxial,  cylindrical  plates  that  carry  current 

in  the  opposite  direction.  The  axes  of  these  lines  coincide  with  the 

major  axis  of  the  toroidal  vessel  and  they  are  located  symmetrically 

around  the  minor  axis  of  the  torus.  These  transmission  lines  change 

mainly  the  local,  vertical  magnetic  field,  while  the  magnetic  flux 

through  the  beam  orbit  remains  approximately  constant.  The  mismatch 

between  the  beam  energy  and  the  vertical  field  is  alleviated  by  a  strong 

focusing  field.  This  field  is  generated  by  a  set  of  t.  ■  2  torsatron 

windings,  i.e.,  two  twisted  wires  that  carry  current  in  the  same 

direction.  In  addition  to  the  transverse  components  of  the  field,  the 

torsatron  windings  provide  a  zero  order  toroidal  magnetic  field.  The 

purpose  of  the  resistive  chamber  wall  is  to  facilitate  the  beam 
2 1 

trapping  in  the  applied  magnetic  field.  Beam  capture  in  the  reba- 
trons  is  very  difficult,  because  the  strong  focusiftg  field  makes  the 
particle  orbit  insensitive  to  the  energy  mismatch  and  thus  small  changes 
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in  Che  betatron  (vertical)  magnetic  field  are  not  sufficient  to  move  the 
beam  from  the  injection  position  near  the  wall  to  the  minor  axis  of  the 
torus. 

The  superior  confining  properties  of  twisted  quadruple  fields  have 
been  recognized  for  several  years.22  Recently,  it  has  been  reported19 
that  the  1*2  Stellarator  configuration  has  an  energy  bandwidth  of 
-50%. 

In  this  report,  we  are  presenting  results  from  our  studies  of  the 
beam  dynamics  in  a  rebatron  accelerator  when  the  magnetic  fields  are  not 
a  function  of  time.  In  addition,  the  local  vertical  magnetic  field  has 
been  replaced  by  a  betatron  magnetic  field.  Our  results  indicate  that 
both  the  bandwidth  and  the  maximum  electron  current  that  can  be  confined 
by  a  rebatron  are  very  high.  Results  with  the  local,  fast  varying 
vertical  field  will  be  reported  in  a  forthcoming  publication. 


II.  The  Applied  Fields 
a.  Magnetic  Fields 

AAA 

In  the  local  cyclindrical  coordinate  system  e^ ,  e^ ,  eg  shown  in 
Fig.  2,  the  magnetic  field  components  of  the  1*2  torsatron  are  given 
by 
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components  in  the  numerical  integration  of  the  orbit  equations.  It  has 
been  determined  that  the  first  two  non  zero  terms  in  the  expansion  are 
sufficient  to  describe  the  field  in  the  region  o/a  <  0.5  with  an 
accuracy  better  than  95%. 

In  the  analytical  work  described  in  Section  IV,  the  toroidal 
corrections  have  been  neglected  as  well  as  all  the  terms  with 
m  >  2.  Furthermore,  it  has  been  assumed  that  *«1.  Under  these 
simplifications  the  torsatron  magnetic  field  becomes 

Bq  =  2B®Xetl2  (x)  sin  [2(A  -  ots)],  (9) 
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For  x  «  l,  Eqs.  (9),  (10)  and  (11)  become 
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In  addition  to  the  torsatron  field,  the  rebatron  accelerator 

includes  a  betatron  or  vertical  magnetic  field  and  a  toroidal 
ex 

field,  B  ,  that  is  produced  by  a  set  of  toroidal  coils. 


The  two 


components  of  the  betatron  field  are  described  by  the  linearized 
equations 
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v.’Vre  B  is  the  betatron  field  at  the  reference  orbit,  i.e.,  at 
zo 

x  =>  y  =  0  and  n  is  the  external  field  index. 

The  toroidal  field  produced  by  a  set  of  toroidal  coils  is 
independent  of  toroidal  angle  and  therefore  varies  as  1/r.  This 
toroidal  field  can  be  chosen  to  have  either  the  same  or  opposite 
polarity  to  the  torsatron  toroidal  field. 


b.  The  Electric  Field  in  the  Gap 

Consider  two  cylinders  with  their  axes  lying  along  the  same  line 
and  separated  by  a  distance  d  as  shown  in  Fig.  4a.  Since  the  cylinder 
on  the  left  is  charged  to  -  VQ  and  the  cylinder  to  the  right  is  charged 
to  +  VQ  the  average  electric  field  in  the  gap  is  <Eg>  ■  2  VQ/d  .  The 
local  electric  field  is  given  by  the  solution  of  Laplace  equation,  i.e., 
«  =  0  .  For  |s|  >  d/2  ,  the  exact  components  of  the  electric  field 

are: 
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Similarly,  for  s  <  d/2  ,  the  two  components  of  the  electric  field 
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where  J  (X  a)  =  0,  a  is  the  radius  of  the  cylinders  and  J„  and  J,  are 
on  o  I 

the  Bessel  functions. 

The  electric  field  lines  that  correspond  to  the  field  components 

given  by  Eqs.  (17)  to  (20)  are  plotted  in  Fig.  4b.  These  electric 

fields  are  a  good  representation  of  the  fields  produced  inside  the  torus 

by  a  transmission  line,  since  in  this  region  the  inductive  magnetic 

2 

field  is  zero  and  therefore  the  potential  is  described  by  17  4>  *  0  . 


III.  Numerical  Results 


To  investigate  the  confining  properties  of  the  torsatron  magnetic 
field,  we  have  integrated  the  relativistic  equations  of  motion  using 
Eqs.  (1)  to  (7)  for  the  torsatron  magnetic  field  and  Eqs.  (15)  and  (16) 
for  the  betatron  field.  The  accelerating  gap  is  2  cm  wide  and  as  shown 
in  Fig.  5,  the  electric  field  is  limited  to  a  0.60  radian  wide  toroidal 
sector.  For  reasons  that  are  discussed  later  on,  the  self  fields  have 


been  omitted  in  these  runs. 


In  the  first  run,  the  current  in  the  torsatron  windings  is  chosen 


to  be  zero.  Figure  6a  shows  the  normalized  particle  energy  (y)  as  a 
function  of  time  and  Fig.  6b  the  projection  of  the  particle  orbit  in  the 
transverse  plane.  The  various  parameters  for  this  run  are  listed  in 
Table  II.  Since  e  is  zero,  the  magnetic  field  configuration  is  that  of 
the  modified  betatron.  As  a  consequence  of  the  curvature  drift,  the 
gyrating  particle  drifts  out  of  the  system  in  about  26  nsec,  i.e.,  in 
about  a  revolution  around  the  major  axis.  As  expected,  the  guiding 
center  of  the  particle  moves  mainly  in  the  vertical  direction,  while  the 
particle  gyrates  around  its  guiding  center  with  a  frequency 
corresponding  to  the  local  toroidal  field. 

Figures  7a  and  7b  show  the  normalized  energy  of  the  particle  and 
its  orbit  when  approximately  -124.7  kA  of  current  flows  through  the 
torsatron  windings.  The  rest  of  the  parameters  for  this  run  are  listed 
in  Table  III.  The  particle  remains  confined  for  eight  revolutions. 
Figures  8a  and  8b  show  similar  results  when  the  current  in  the  windings 
is  increased  to  approximately  -250  kA.  The  corresponding  torsatron 
field  strength  factor  e  is  -0.8.  The  remainder  of  the  parameters  are 
listed  in  Table  IV.  In  all  three  runs  the  betatron  magnetic  field  was 
held  constant  at  118  G.  These  results  clearly  demonstrate  that  the 
confining  properties  of  the  system  are  substantially  improved  by  the 
addition  of  the  torsatron  field.  The  particle  strikes  the  chamber  wall 
when  its  gamma  approaches  approximately  65.  The  total  time  the  particle 
remains  in  the  system  is  about  320  nsec,  i.e.,  more  than  an  order  of 
magnitude  longer  than  when  the  torsatron  field  is  absent. 

Further  improvement  in  the  particle  confinement  is  observed  when 
the  period  of  the  windings  is  reduced  or  the  current  in  the  windings 
increased.  An  additional  modest  improvement  in  the  confinement  of  the 
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system  Is  observed  when  the  betatron  field  is  increased  above  its 
matching  value.  This  is  shown  in  Fig.  9.  The  betatron  field  for  this 
run  is  236  G  and  the  remainder  of  the  parameters  are  identical  to  those 
in  Fig.  8.  The  confinement  time  increased  by  20  nsec,  i.e.,  from  320  to 
340  nsec.  However,  when  the  betatron  field  increased  to  472  G  the 
confinement  time  was  reduced  to  290  nsec. 

IV.  Theoretical  Model 

To  gain  a  better  understanding  of  the  focusing  properties  of  the 
torsatron  fields,  we  have  developed  a  theoretical  model  that  is  based  on 
linear  external  fields.  Obviously,  these  fields  are  appropriate  only 
near  the  minor  axis  of  the  torus,  i.e.,  when  2oto  «  1. 

The  components  of  the  torsatron  field  in  the  coordinate  system 
er’  e4’  ez  shown  in  FiS*  2  are 

B  =  B  cos -  B,  sinii  ■  BeXe«  ^zcos2nr  4  +  (r-r  )  sin  2ar  4 1  (21) 

rt  o  st  o  o  o 
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In  addition,  the  betatron  magnetic  field  is  given  by 
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where  n  is  the  external  field  index.  The  total  field  components  are 
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GX 

where  Bg  indicates  any  additional  toroidal  field  that  may  be  applied. 
The  accelerating  electric  field  components  are  approximated  by 
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where  y  =*  dy/dt  ,  y  *  d  y/dt  and  v  is  the  toroidal  velocity,  which  is 
assumed  constant. 

Using  Eqs.  (26)  to  (31)  for  the  fields,  the  equations  of  motion  in 
the  laboratory  frame  become 
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Since  at  the  start  of  the  acceleration  oi  oi  /2  »  *ui  and 
0  low  l  I  l 

(y/2y)  |( — - — )|  »  |  Soi2  |  ,  the  two  coupled  equations  (36)  and  (37)  can 
be  combined  into  a  single  equation  by  introducing  the  complex  variable 
’!i  =  R  +  i  Z  .  Multiplying  Eq.  (37)  by  i  and  adding  it  to  Eq.  (36),  it 
is  obtained 
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Equations  (36)  and  (37)  have  been  solved  numerically.  After  integration 
the  orbit  is  transferred  back  to  the  laboratory  frame.  The  results  are 
shown  in  Fig.  10.  The  projection  of  the  orbit  in  the  r,z  plane  is  shown 
in  Fig.  10a,  the  particle  radial  distance  from  the  minor  axis  as  a 
function  of  time  is  shown  in  Fig.  10b  and  y  as  a  function  of  time  in 
Fig.  10c.  The  various  parameters  for  this  run  are  identical  to  those 

listed  in  Table  IV.  The  particle  strikes  the  wall  at  about  325  nsec, 

when  its  gamma  is  approximately  68.  These  results  are  in  good  agreement 
with  those  of  Fig.  8  that  have  been  obtained  using  the  more  accurate 

expressions  for  the  torsatron  fields.  As  will  be  discussed  later,  the 

particle  was  lost  because  at  Y  a  65  it  entered  the  unstable  region  that 

extends  from  y  =  65  to  y  =  121  . 

When  y  =  0  ,  the  homogeneous  part  of  Eq.  (38)  becomes 

*  +  [2fj  +  i\)  Hi  +  [fj  -  f*]  *  -  0,  (39) 

i.e.,  a  fourth  order  equation  with  constant  coefficients.  The  solutions 

of  Eq.  (39)  are  of  the  type  '!>  *  e^  Jt  ,  where  ,  is  given  by 


2 

when  (d;  «  a)  a)  /2 

i  o  w 


The  particle  orbits  are  stable  when  <d  is  real,  i.e.,  when 


u(y)  s 


>  0. 


(41) 


The  two  roots  of  (40)  are  given  by 


Y 


± 


2% 


U 

w 


f  1  ± 


(C 

t 


n 


1  . 


(42) 


The  function  u(y)  is  plotted  in  Fig.  11  for  three  values  of  B  and 

n 

B  6X.  In  Fig.  11a  B  ex  -  -6kG,  B  -  5KG  and  thus  BQ  *  -(B  eX  +  Brt)  =  1KG. 
s  so  8  s  0 

For  this  value  of  toroidal  magnetic  field  and  for  »  0.4  and 

9  -1 

o>w/2  =  3x10  sec  ,  Eq.  (42)  gives  19.96  and  -8.2  for  the  two  roots  of 
Eq.  (41).  Therefore,  at  y  ■  7  the  particle  orbit  should  be  unstable. 
Results  from  the  numerical  integration  of  nonlinear  orbit  equations  for 
Y  =  0  ,  and  Y  *  7  and  using  the  same  values  for  the  rest  of  the 
parameters  as  in  Fig  11a  are  shown  in  Fig.  12.  As  expected,  the  orbit 
is  indeed  unstable  and  the  particle  is  lost  in  less  than  one  nsec. 

By  reversing  the  direction  of  the  current  in  the  torsatron  wires 
Bq  and  et  change  sign  and  the  two  roots  of  Eq.  (41)  become  50.5  and 
78.7.  Therefore,  for  y  *  7  the  orbit  is  stable.  This  is  in  agreement 
with  the  results  from  the  numerical  integration  of  nonlinear  orbit 
equations  shown  in  Fig.  13. 

When  y+  <1  ,  the  orbits  are  stable  for  all  values  of  Y  •  For 
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6X 

B  *  6KG  ,  B  =  5KG  and  e  *  -0.4  the  two  roots  of  Eq.  (41)  are 
sot  ^ 

-50.5  and  -78.7.  For  this  case  the  orbits  were  found  stable  for  all  the 
values  of  v  considered. 

The  numerical  and  theoretical  results  are  in  excellent  agreement 

0X 

when  the  linearized  theoretical  model  is  valid.  However,  when  B  and 

s 

Bq  have  opposite  signs  and  |BgeX|  >  |B0|  •  the  toroidal  field  vanishes 
at  some  radial  distance  and  the  field  lines  form  magnetic  cusps.  In 
this  case  the  linear  theory  does  not  properly  describe  the  fields  and 
the  predictions  of  the  theory  are  not  in  agreement  with  the  numerical 
results. 


When 


r"o  “» /(»w  -r)2  <\  >’  «  *  • 

Y 


Equation  (40)  gives 


(m  nr  ri  + 

w 


VY>(Vy) 


1 ,  (slow  mode) 


“+  = 


w 


2  2 

'J)  u 

o  w 


w  Y 


4(u>w  -  H^/yXO^/y) 


(fast  mode) 
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In  the  laboratory  frame  the  slow  mode  0.  becomes 


The  particular  solution 

of  Eq.  (38) 

in  the  rotating  frame 

2 

iV  V 

=  0  and  small,  is 

V  =*  R  + 

iZ 

P  P 

P 

where 

4A2(oi  -  0  /y 

+  'oi  /  2 ) 

(0  t 

~  T*7  A 

R  =*  - — - 

0 

cna( — — - 1 

Transforming  back  to  the  laboratory  frame  using  the  transformation 

-i  'o  t/2 
~  w 

®  “  e  ,  we  find  that  the  particle  orbit  is  displaced  along  the 

horizontal  axis  by 


to  y 
o 


%/y,V 


Figure  14  shows  the  projection  of  the  particle  orbit  in  the 
tranverse  plane  for 

9  -1 

Y  *  11,  e  -  -0.4  B.  =*  11KG  ,  »o  *  6  x  10  sec  ,  r  *  100cm, 
t  °i  w  ’  o 

B  n  *  118  G  and  u.  »  c.  For  these  parameters  Eq.  (43)  gives  a  slow 


period  t+  =■  2ir/°+  -  62  nsec.  For  the  same  parameters  the  code  gives 

t+  =  60nsec. 

In  addition,  Eq.  (46)  gives  a  displacement  Ar  =  1.74  cm  ,  which  is 

identical  to  the  orbit  displacement  of  Fig.  14. 

Let’s  now  return  to  discuss  briefly  the  results  of  Fig.  10.  For 

the  parameters  of  the  run,  Eq.  (42)  gives  y+  3  121  and  Y_  *  65  .  When 

the  y  of  the  particle  reaches  65  i.e.,  at  about  300  nsec,  it  becomes 

unstable  and  strikes  the  wall  in  one  revolution. 

In  addition,  at  t«0  the  ratio  0 q/oj  y  =  6.67  and  according  to  Eq. 

(46)  the  orbit  displacement  is  negative.  As  y  increases  Q  /y m  is 

B  w 

reduced  and  when  /yu  <  1  the  orbit  dispaceraent  becomes  positive.  At 

“  w 

y  =  46.6  ,  0  /y'o  3  1  and  Ar  =  0  .  According  to  Fig.  10b  this  occurs 

H  W 

at  t  3  210  nsec  ,  which  corresponds  to  y  3  47  (see  Fig.  10c). 

V.  Self  Fields 

An  accurate  self  consistent  determination  of  self  fields  of  a  high 
current  electron  ring  confined  in  a  rebatron  is  difficult,  because  the 
minor  cross  section  of  the  ring  has,  in  general,  a  complex  shape  that 
varies  along  the  toroidal  direction. 

Since  we  are  interested  in  the  macroscopic  motion  of  the  ring  and 
therefore  on  the  self  fields  that  act  on  the  ring  centroid,  we  assume 
that  the  ring  has  a  circular  cross  section  and  its  particle  density  is 
uniform.  Neglecting  toroidal  corrections,  the  fields  at  the  center  of 

the  beam,  which  is  located  at  the  distance  (r  -  r  )  and  z  from  the 
18,26 


minor  axis  are 


E 

r 


(47) 


*  2*lel  Vo  “T 

a 


-2  w  |  e  |  n  r 

q  o  L  r 
a  o 


(48) 


BS  *  2ir|e|n  8  r 
r  !  I  o  o  i 


_b  _z 

o  2  r 
a  o 


(49) 


and 


K  <r  -  V 

2ir  I  e  I  n  «  r  — -  - 

I  '  o  o  o  2  r 


(50) 


where  nQ  is  the  particle  density,  r^  the  beam  radius  and  a  the  minor 


radius  rf  the  perfectly  conducting  torus. 

When  y  *  0,  n  =  V2  and  the  beam  energy  is  matched  to  the  vertical 
field,  the  equation  describing  the  beam  centroid  motion  in  the  tranverse 
rotating  plane  is  given  by 


"  -  ^  _  *  ~  2 
'!)  +  f  j  4  +  f  —  if^’^3*  Xe 


(51) 


where  f^  and  f^  have  been  defined  under  Eq.  (38)  and 


:  -2  ,2  .  “»  % 


(52) 


where 


*2 


2  r  2y 


i_  (ib,2 

3  ka  ;  » 


(53) 


and  'i>.  ■  4ire  n  /m  . 


iu) .  t 


2  ~  _ 

The  solution  of  Eq.  (51)  when  X  *  0  is  'b  =  v  e  , 


where 


(-T  +  * 


4y 


o  2  2  h 

+  r  ,  QS  v  2  A2  “o  ^ 

±r  c-  -— )  (-J  +  «L  )  +  — 

4y 


o  w 


(54) 


The  orbits  are  stable  provided 


q  ~2 

—T~  +  hi  >0  , 


(55) 


4y 


and 


Cl) 


03 


r;  2-  (-^-)2  +  -2  —  - 

I  2  2  y 


2  2 

(0  'ii 


’Art  w  xj 

w  a  ,2  o  w 


>0 


Equation  (56)  can  be  written  as 


4  .  A2 
L  +  2“l 

and  its  roots  are  given  by 


~4  .  "2  /  Six,  4  w  x  .  ,  w  >2r,  b  w.2  z  i  . 

u  -  uL  +  2ml  (—)(“ - j  5  - 2  “  %  >0  » 


A2 

0), 


9.  w 

(»/»-- r<r - r>  *V 


w 


(56) 


(57) 


(58) 


°o  01 

When  |( - j)  |  >  |  u>o  |  ,  the  two  roots  of  inequality  (57)  are  both 


either  positive  or  negative  depending  upon  the  sign  of 


•">  fjj 

^  w 

( - — )  and  a)  .  The  results  are  summarized  in  Table  V. 

y  2  o 


For  y»l,  Eq.  (58)  becomes 


i.e.  the  orbits  are  always  stable  provided  inequality  (55)  is  satisfied. 

The  maximum  amount  of  space  charge  that  can  be  confined  by  a 
rebatron  can  be  determined  from  Eq.  (55).  When 


« 

,,  5  a 

2y 


Eq.  (55)  gives 


_  2  2 
%  a 


8c 


(60) 


For  a  =  10  cm,  B.=>  10  KG,  y  =*  7,  Eq.  (60)  gives  v  °>  3,000  or  I  *  50  MA. 


"2 

When  the  current  of  the  beam  is  large,  «0.  However,  as 

*2  2  2 
y  increases  w.  approaches  it  asymptotic  value  u  /2r  .  Similarly,  the 

"2 

two  roots  m|  +  approach  their  asymptotic  value  given  by  Eq.  (59). 


id 

Figure  15  shows  the  stability  diagram  at  t  *  0,  when  - ^  0 

and  o>o  <  0.  For  this  case  the  stability  condition 


is  - 


-J£  r(_i 

2  Y 


1 1) 

-2)1 

2 


'0  (x) 

O  w  .  "2 
+  — ; —  >  0) 

2  I 


>  - 


'Pi 


4y 


2  2  ^  2 

During  acceleration  y  increases  and  therefore  both  -°.  / 4y  and  «j  move 

to  the  right  of  the  diagram.  Therefore,  it  may  be  argued  that 

*  2  2  2 

before  crosses  the  vertical  axis,  /4y  catches  up  with  it  and  the 

"2 

ring  becomes  unstable.  A  similar  situation  would  occurs  when  w  becomes 
A2 

equal  to  u>|  +  .  However,  we  have  shown  that  when  the  system  is  stable  at 
t  *  0  it  will  remain  stable  for  any  y  that  exceeds  the  initial  y  . 
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VI.  Conclusions 


We  have  carried  out  an  extensive  numerical  and  analytical 
investigation  of  the  beam  dynamics  in  a  rebatron  accelerator.  Although 
the  analytical  work  is  based  on  simple,  linear  approximations  for  the 
various  fields,  the  two  approaches  give  very  similar  results  when  these 
approximations  are  valid. 

Our  studies  indicate  that  when  self  field  effects  can  be  ignored, 
the  particle  normalized  energy  can  be  increased  from  y  =  7  to  y  =*70,  at 
constant  betatron  field,  before  confinement  is  lost.  This  implies  that 
the  device  has  a  bandwidth  that  approaches  1000%.  This  bandwidth  can  be 
further  increased  by  increasing  the  current  in  the  torsatron  wires. 

Even  in  the  absence  of  the  space  charge,  there  is  a  range  of 
parameters  [see  Eq.  (42)]  for  which  the  rebatron  is  unstable.  However, 
this  orbit  instability  can  be  easily  avoided  by  a  judicious  choice  of 
the  various  parameters. 

As  far  as  orbit  stability  is  concerned,  the  maximum  electron  beam 
current  that  can  be  confined  in  a  rebatron  accelerator  is  given  by  Eq. 
(60)  and  is  impressively  high.  Therefore,  it  is  expected  that  the 
limiting  beam  current  in  a  rebatron  would  be  determined  from  collective 
instabilities  and  not  from  the  macroscopic  stability  of  beam  orbits. 

Although  the  bandwidth  of  rebatron  accelerators  is  very  high,  the 
maximum  energy  that  can  be  obtained  by  these  devices,  with  time 
independent  magnetic  fields,  is  rather  limited.  To  achieve  very  high 
energies  (y  >  1000)  the  betatron  magnetic  field  should  be  replaced  by  a 
local  vertical  magnetic  field  that  varies  rapidly  with  time  and 
approximately  in  synchronism  with  the  beam  energy.  Such  a  fast  vertical 
field  can  be  generated  by  two  coaxial,  cylindrical  lines  that  carry 
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current  in  the  opposite  direction.  The  axes  of  these  lines  coincide 
with  the  major  axis  of  the  toroidal  vessel  and  they  are  located 
symmetrically  around  the  minor  axis  of  the  torus.  These  transmission 
lines  change  mainly  the  local  vertical  magnetic  field,  while  the 
magnetic  flux  through  the  beam  orbit  remains  approximately  constant. 
The  mismatch  between  the  beam  energy  and  the  vertical  field  is 
alleviated  by  the  strong  focusing  field.  The  effect  of  the  rapidly 
varying  vertical  magnetic  field  on  the  beam  dynamics  will  be  reported  in 


a  forthcoming  publication 


Table  I 


Parameters  relevant  to  the  torsatron  fields  shown  in  Fig.  3.  Only  two 
terms  retained  in  the  series  of  Eqs.  (1)  to  (3). 


Torus  major  radius  r  (cm)  =  100 

o 

Windings  minor  radius  0Q  (cm)  =  12 

Toroidal  chamber  minor  radius  a  (cm)  =  10 

a  =  2-^/L  (cm  1 )  =0.1 

Field  Strength  Factor  =  0.2 

Winding  Current  I(kA)  =  62.37 

e  =  2 

ex 

Additional  Toroidal  Field  B  (kG)  =  -  6 


Table  II 


Parameters  of  the  run  shown  in  Fig.  6 


Torus  Major  Radius  r0(cm) 

Winding  Minor  Radius  oq  (cm) 
Toroidal  Chamber  Minor  Radius  a(cm) 
a  =  2ir/L  (cm  *) 

Field  Strength  Factor  s 
Winding  Current  I  (kA) 

1 

gy 

Additional  Toroidal  Field  Bg  (kG) 

Betatron  Field  B  (  G) 

zo 

Ext.  Field  Index  n 


Initial 

Initial  Positions 
Initial  Velocities 


Y 


0 

V 


<*)  =  9 


Table  III 


Parameters  of  the  run  shown  in  Fig.  7. 


Torus  Major  Radius  rQ(cm) 

Winding  Minor  Radius  pQ  (cm) 

Toroidal  Chamber  Minor  Radius  a(cm) 
a  =  2tt/L  (cm  * ) 

Field  Strength  Factor  e 
Winding  Current  I  (kA) 

?. 

Additional  Toroidal  Field  BeX  (kG) 

s 

Betatron  Field  Bzq(G) 

Ext.  Field  Index  n 

Initial  y 

Initial  Positions  o  *  <b  ■  s 

Initial  Velocities  v  =  v  =  0  ,  v 


Table  IV 


Parameters  of  the  run  shown  in  Fig.  8. 


Torus  Major  Radius  rQ(cm)  =  100 

Winding  Minor  Radius  po(cm)  =  12 

Torodial  Chamber  Minor  Radius  a( cm)  *  10 

a  =  2tt/L  (cm  *  )  =0.1 

Field  Strength  Factor  e  ■  -0.8 

Winding  Current  I  (kA)  *  -250 

l  =2 

Additional  Torodial  Field  BSX(kG)  =  -6 

s 

Betatron  Field  B  (G)  =  118 

zo 

Ext.  Field  Index  n  =  0.5 

Initial  y  =7.0 

Initial  Positions  o  =  <b  =  s  =0 

Initial  Velocities  v  **  v  =  0  ,  v  »  c 
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a 


both  negative 


both  positive 


both  positive 


both  negative 


lematic  of  a  rebatron.  The  accelerating  gap  is  energized  by  two  or  more  transmission  lines 
that  are  symmetrically  distributed  around  the  minor  axis. 


Fig.  3  —  Torsatron  magnetic  field  components  (a)  at  </>  *■  s  “  0  and  (b)  at  s  —  0,  <t>  ■■  n/2.  In 
addition  to  the  torsatron  field  there  is  toroidal  field  5“  “  —6  kG  that  is  produced  by  a  set  of  toroidal 
coils,  (c)  magnetic  field  lines  in  r,  z  and  p,  s  planes. 


Fig.  3  (Cont’d)  —  Torsatron  magnetic  field  components  (a)  at  <S>  *■  s  —  0  and  (b)  at  s  ™  0,  <f>  —  n/2. 
In  addition  to  the  torsatron  field  there  is  toroidal  field  Bfx  —  -6  kG  that  is  produced  by  a  set  of 
toroidal  coils,  (c)  magnetic  field  lines  in  r,  z  and  p,  s  planes. 


Fig.  3  (Cont’d)  —  Torsatron  magnetic  field  components  (a)  at  4>  —  s  *■  0  and  (b)  at  s  —  0,  <t>  «■  n/2. 
In  addition  to  the  torsatron  field  there  is  toroidal  field  5“  —  -6  kG  that  is  produced  by  a  set  of 
toroidal  coils,  (c)  magnetic  field  lines  in  r,  z  and  p,  s  planes. 
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(b) 

Fig.  10  —  (a)  Particle  orbit  in  the  r,  z  plane;  (b)  particle  radius  as  a  function  of  time  and  (c)  y  of 
particle  as  a  function  of  time.  These  results  have  been  obtained  from  the  linear  equations  (36)  and 
(37).  The  results  shown  are  in  the  Lab.  frame.  The  various  parameters  for  this  run  are  the  same  with 
those  of  Figure  8. 


Time  (nsec) 
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Fig.  10  (Cont’d)  —  (a)  Particle  orbit  in  the  r,  z  plane;  (b)  particle  radius  as  a  function  of  time  and  (c) 
y  of  particle  as  a  function  of  time.  These  results  have  been  obtained  from  the  linear  equations  (36)  and 
(37).  The  results  shown  are  in  the  Lab.  frame.  The  various  parameters  for  this  run  are  the  same  with 
those  of  Figure  8. 


Particle  Orbit  in  The  Transverse  Plane 
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Fig.  14  —  Particle  orbit  in  the  r,  z  plane  for  an  initially  mismatched  beam  and  y  -  0.  The  particle 
rotated  around  the  equilibrium  position  four  times  with  a  period  of  60  nsec.  The  theory  [Equation 
(43)]  predicts  a  period  of  62  nsec  and  the  nonlinear  equations  give  similar  orbits  with  a  period  of  about 
55  nsec. 
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